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The evolution of scalar fields transported by turbulent flow is characterized by the presence of fronts, 
which rule the small-scale statistics of scalar fluctuations. With the aid of numerical simulations, 
it is shown that : isotropy is not recovered, in the classical sense, at small scales; scaling exponents 
are universal with respect to the scalar injection mechanisms ; high-order exponents saturate to a 
constant value ; non-mature fronts dominate the statistics of intense fluctuations. Results on the 
statistics inside the "plateaux" , where fluctuations are weak, are also presented. Finally, we analyze 
the statistics of scalar dissipation and scalar fluxes. 



I. INTRODUCTION 

The understanding of small-scale fluctuations in scalar 
fields, such as temperature, pollutant density, chemical 
or biological species concentration, advected by turbulent 
flow is of great interest in both theoretical and practical 
domains |l| . Propagation phenomena of light beams and 
radio waves in the atmosphere are for example strongly 
influenced both by the magnitude and the spatial dis- 
tribution of small-scale temperature gradients |g] . Their 
dynamical properties have been a subject of very accu- 
rate experimental investigations carried out in the last 



in the ocean 

A striking feature 



few years both in the atmosphere J 2 
and for laboratory turbulent flow |6 
of all these situations is the presence of fronts (also called 
sheets or cliffs). The important point is that the scalar 
field has very strong variations across the fronts, sep- 
arated by large regions ("ramps" or "plateaux") where 
scalar fluctuations are weak (see for example Fig. pj\ ). 

In the atmosphere, "ramp-and-cliff" structures in the 
temperature field extend from the ground to the middle 
stratosphere M. They are of major importance for the 
ubiquitous stratospheric radar echoes |10|| and the prob- 
lem of their formation is known as frontogcnesis pi| . 
For the temperature field it is crucial to take into ac- 
count its active nature, i.e. the fact that it influences the 
advecting velocity fleld. Simplifled quasi and/or semi- 
geostrophic models have been proposed and investigated 
for the frontogcnesis and we refer the interested reader 
to Refs. [0,Q and references therein. 

Cliffs such as those observed in the active case have in 
fact been observed also for passive situations, i.e. when 
the scalar field does not affect the advecting velocity 
field. Remarkably, cliffs are not mere footprints of ve- 
locity structures and they arise both for velocity fields 
that are solutions of the Navier-Stokes (NS) equations 
p3| , p^ and for synthetic velocity fields [p^-p^ , including 
when the velocit y is rapidly varying in time (the so-called 
Kraichnan case |19|). As emphasized in Refs. ||l3| , |l5| , p0| , 
the velocity gradient stretching is weak in the elliptic re- 
gions of the velocity field, where their 'rotational' char- 



acter inhibits the formation of strong scalar gradients. 
Conversely, compressions along specific directions take 
place in the hyperbolic regions of the incompressible flow, 
scalar particles coming from distant regions of space ap- 
proach each other and strong scalar gradients are thus 
developed. 

The observation of ramp-and-cliff structures naturally 
raises the issue of their consequences for the scalar statis- 
tical properties and their possible implications for mod- 
els of scalar transport. Two points will be specifically 
addressed here for the passive scalar case. Part of the 
results were previously summarized in Ref. |14| . 

The first concerns the role of ramps and cliffs for 
the scalar anisotropy. A central postulate of the clas- 
sical Kolmogorov-Obukhov-Corrsin theory (see, e.g., 
Ref. [2^) for passive scalar turbulence is that the 
anisotropy degree should decay as smaller and smaller 
scales are considered. When the system is driven by an 
anisotropic forcing, such as a large-scale scalar gradient, 
the most common experimental situation, isotropy is usu- 
ally supposed to be fully recovered at small scales. How- 
ever, the evidence stemming from the experiments is that 
this is not the case (at least not in the original sense, see 
the following) . Ramps and cliffs are preferentially aligned 
with the direction of the large-scale scalar gradient and 
this permits the anisotropics to find their way down to 
the smallest excited scales. In particular, the probability 
density function (pdf) for the scalar derivatives along the 
direction of the large-scale scalar gradient are strongly 
skewed. Furthermore, the scalar gradient skewness (that 
should vanish for an isotropic field) is 0(1), indepen- 
dently of the Peclet number ||-||,|l^,|l^ . This leads us to 
investigate scalar turbulence universality, i.e. its degree 
of independence of injection mechanisms. 

Recent numerical and experimental works exploit- 
ing symmetries under rotations point to the fact that 
anisotropic turbulent fluctuations are subdominant with 
respect to isotropic ones. We shall discuss the relation 
between these results and the aforementioned experi- 
mental observations about the persistence of large-scale 
anisotropics. 



The second issue associated to the presence of chffs 
concerns their role for the statistics of scalar differ- 
ence strong fluctuations, i.e. to the behavior of high- 
order structure functions. In particular, scalar variations 
across the cliffs are comparable to the root-mean-square 
(rms) value of the scalar. This suggests that scalar struc- 
ture function scaling exponents might saturate, that is 
tend to a constant for large enough orders. The ques- 
tion of saturation was first raised in Ref. iQ. Here, we 
shall provide evidence for saturation in two different sit- 
uations : scalar advection by a two-dimensional Navier- 
Stokes velocity field in the inverse cascade regime and 
in the three-dimensional Kraichnan model [|l9|. As dis- 
cussed in the sequel, the latter model represents the least 
favorable case to observe saturation, because of the ve- 
locity short correlation time. The fact that saturation is 
still observed points to the genericity of the phenomenon 
for scalar turbulence. 

The paper is organized as follows. In Sections || 



and III, we recall some well-known results about the 



expected classical laws for the decay of anisotropies 
in scalar turbulence and the properties of the Navier- 
Stokes two-dimensional inverse energy cascade, respec- 
tively. Anisotropies and scalar turbulence universality 
are discussed in Section IV, while the results about sat- 
uration are presented in Section [V|. The successive Sec- 
tion investigates ramp- and- cliff structure dynamical pro- 
cesses. The aim is to clarify whether the observed satura- 
tion arises from 'mature' cliffs, having thicknesses compa- 
rable to the dissipative scale, or from 'non-mature' ones, 
still in the process of steepening. Section VII reports re- 
sults on the scalar statistics inside t he "pl ateaux" . Scalar 
dissipation is the subject of Section |VIII , and the results 
on scalar fluxes, involving joint velocity-scalar correla- 
tions, are discussed in Section 
devoted to conclusions. 



IX. The last Section is 



The correlation x {''"/Lf) is concentrated at the forcing 
integral scale Lf and rapidly decreases for r ^ Lj. 
An alternative injection mechanism, closer to typical 
experimental situations, is to maintain a mean scalar 
gradient g. Scalar fluctuations 9 superimposed on it, 
T{r, t) = 9{r, t) + g ■ r, obey the equation of motion : 

dtOir, t) + v{r, t) ■ V e{r, t) = KA0{r, t)-g-v. (3) 

Note that the mean gradient injection mechanism selects 
by its very definition a preferential direction in the system 
and isotropy is thus broken. 

We shall be interested in the passive scalar structure 
functions 

5„(r) ^ {[T{x + r, t) - T{x, t)]") ^ ((-^.Tf), (4) 

where ( ) denotes the average over the velocity, and pos- 
sibly the forcing, ensemble. The reason to focus on cor- 
relations of the field T rather than those of the field 6 
is explained in Appendix Kl Throughout this work we 
assume the velocity to be scale-invariant, at least in the 
scalar inertial range of scales defined by ry <C |r| <C L. 
Here, 77 is the scalar dissipation scale and L is the scalar 
integral scale. For the mean-gradient injection, L is com- 
parable to the velocity integral scale. For the randomly 
forced case, we shall consider situations where the veloc- 
ity integral scale is comparable (or even larger) than the 
forcing one and therefore L ~ Lf. On account of these 
assumptions, the scalar structure functions are expected 
to have a scaling behavior in the inertial range of scales. 
For the isotropic randomly forced case, odd orders vanish 
by symmetry. For the injection by a mean gradient, both 
even and odd-order structure functions are non-zero. 

For even-order structure functions, the classical predic- 
tion of the Kolmogorov-Obukhov-Corrsin (KOC) theory 



II. THE CLASSICAL THEORY OF SCALAR 
TURBULENCE 

The advection-diffusion partial differential equation 
governing the evolution of a passive scalar field T{r,t) 
is: 



dtTir, t) + v{r,t)-V T{r, t) = KAr(r, t). 



(1) 



Here, t)(r, t) is the incompressible advecting velocity field 
and K is the molecular diffusivity. The "energy" T^ is 
statistically conserved by the advection term in ([I]) and 
dissipated by the viscous term. In order to attain a sta- 
tionary state it is thus necessary to inject scalar fluctua- 
tions. The simplest way is to add a forcing term f{r,t) 
to the right-hand side (rhs) of (nl). A convenient choice 
that we shall use in this work is to take the forcing / 
random, Gaussian, statistically homogeneous, isotropic, 
white in time, of zero mean and correlation function 



if {r,t)f (0,0)) =x{r/Lf)6{t). 



(2) 



^2„(r) = C2„6^-"/V"/M- ocr';- 



(5) 



with (^2n = 2n/3 and the anomaly p2n = 2n/3 — (^2n = 0. 
In the previous formula, ee = k((V0)^) is the mean scalar 
energy dissipation, C'2n are non-dimensional constants 
and the velocity is assumed to be of Kolmogorov-type, 
with velocity increments (((5rw)") ex r"/"^ and a finite en- 
ergy flux ey. The same arguments can be reformulated 
along the same lines for other types of velocity fields. 

Dimensional arguments are easily extended to odd- 
order moments in the presence of a mean gradient g P2| . 
The balance between the advection term and the injec- 
tion term on the rhs of (0) gives indeed 
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(6) 



with C2n-i-i = 2n/3 -I- 1. The classical prediction for the 
decay rate of hyperskewnesses Rn is then : 



R„ 



8211 



+1 



"52 



-1/2 



p2/3 



1,2,-. 



(7) 



Predictions for the scaling of scalar gradient moments 
with the Peclet number are obtained by substituting r = 
noc Pe"'^' * in the structure function expressions (^ and 
([7|). It follows in particular that, in the presence of a 
large-scale gradient g, the skewness of the scalar gradient 

1/9 

along this direction should decrease as Pe ' . 

Notice that the classical KOC scalings do not explicitly 
involve the integral scale L. Anomalous scalings corre- 
spond to the violations of these behaviors, i.e. C„ 7^ n/i, 
and the appearance of the integral scale L in the structure 
function expressions via the non- vanishing anomalies pn- 



III. THE ADVECTING VELOCITY FIELDS 

We discuss now the choice of the advecting velocity 
field in Eqs. (0) or (|). The flow that we shall con- 
sider mostly in this work is the one arising from a two- 
dimensional Navier-Stokes inverse energy cascade [ p3[ . 
Its properties are briefly recalled in the following subsec- 
tion. In order to corroborate the results on saturation, 



we shall also consider in Section VB three-dimensional 
synthetic flows of the Kraichnan type (see Refs. |13,E4|). 



A. The 2D Navier-Stokes flow 

The two-dimensional Navier-Stokes equation for the 
vorticity Lo{r^ t) — — A-0(r, t) is : 



dtuj + J (w, ip) = vlS.u} — auj — AT, 



(8) 



where the velocity v is related to the stream function tp 
as tJ = V^V-' — iy yi' ^ ~'^ xi') , the Jacobian is denoted 
by J and the forcing by T. The friction linear term 
—OLijj extracts energy from the system at scales compa- 



1/2 
fit, a 



3/2 



assummg a 



rable to the friction scale 77^. 

Kolmogorov scaling law for the velocity. As before, £„ 

indicates the mean energy dissipation. 

Let us first recall from Ref. pSl some technical details 
on the numerical integration of (8|). To focus on the in- 
verse cascade phase and avoid the Bose-Einstein conden- 
sation in the gravest mode |23|, we choose a to make r^^ 
sufficiently smaller than the box size. The other relevant 
length in the problem is the small-scale forcing correla- 
tion length Ij^, bounding the inertial range for the inverse 
cascade as Z;r <c r ^ r^fr . We use a Gaussian forcing with 
correlation fimction {T{r,t) T{0,t')) = S{t-t') F{r/ljr). 
The forcing correlation decays rapidly for r ^ lyr and 
we choose F{x) = FoZ^exp(— a;^/2), where Fq is the en- 
ergy injection rate. The integration is performed by a 
standard 2/3-dealiased pseudospectral method on a dou- 
bly periodic square domain 2048 x 2048. The viscous 
term in Eq. (@) has the role of removing the enstrophy 
at scales smaller than Ijr and, as customary, it is numer- 
ically more convenient to substitute it by a hyperviscous 



term (of order eight in our simulations). Time evolu- 
tion is implemented by a standard second-order Adams- 
Bashforth scheme. The integration is carried out for one 
hundred eddy turn-over times after both the velocity and 
the scalar fields have reached the stationary state. 

The inverse nature of the energy cascade clearly 
emerges from the third-order longitudinal structure func- 
tion Vii '{r) = {[{v{r,t) — v{0,t)) ■ r/r] ). In analogy 
with the 4/5 law valid for the three-dimensional case (see, 
e.g., Ref. P6|), one can indeed easily derive from (||) the 
relation V} (r) = 3/2evr, expressing the upward con- 
stant energy flux. The neat plateau at the positive value 
3/2 visible in Fig. Ill is a clear confirmation of the inverse 
energy cascade. 




FIG. 1. Compensated third order longitudinal structure 
function V^, (r)/(ei,r). The dotted line is the value 3/2. The 
vertical scale is linear. 

An important property of the flow that has emerged 
both from numerical simulations [ p7| , p5[ and experiments 
||28| is the absence of intermittency in the velocity incre- 
ment statistics. No corrections to the Kolmogorov scal- 
ing r"''^ for the n-th order velocity structure function 
can be measured. In Fig. we show indeed that ve- 
locity increment pdf's at different separations r can all 
be collapsed rescaling the increment S^-v by their vari- 
ance {{Srvyy^"^ '^ r^^^. Furthermore, deviations from 
a Gaussian are small (at least for fluctuations not very 
large) and this has been exploited to propose a theory 
for the observed absence of intermittency |29|. 

The important point to our purposes here is that the 
velocity field arising from the 2D inverse energy cascade 
is scale-invariant (no intermittency) of exponent 1/3, 
isotropic and, mostly important, has realistic correlation 
times. Contrary to synthetic flows, like thos e in psf or 
the Kraichnan ensemble used in Section VB, the corre- 



lation times are finite and free of any sweeping-related 
pathology. The anisotropic and intermittency properties 
of the scalar field that will be discussed in the sequel 
are therefore intrinsic of the advection-diffusion passive 
scalar dynamics and not mere footprints of the advecting 
velocity. 
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FIG. 2. Pdf's of velocity fluctuations rescaled by their 
variance, for five lengthscales within the inverse energy cas- 
cade range. The dashed line is the Gaussian probability den- 
sity. 



IV. ANISOTROPY AND UNIVERSALITY 

The classical picture of the anisotropy decay presented 
in Section O is contradicted by the experimental obser- 
vations. Contrary to the expected Pe~ ' behavior, the 
skewness of the scalar derivative along the large-scale gra- 
dient 3||0 = g ■ Vfl/lgl does not show any measurable 
dependence on the Peclet number ||^-|9|,|30[| . The corre- 
sponding pdf P(9||6') is strongly asymmetric with a sharp 
maximum located at —5, signaling the mean gradient ex- 
pulsion out of the plateaux regions. The same general 
picture is also valid for passive scalar fields advected by 
synthetic flows. 

The fact that the memory of the large-scale injection 
conditions is not lost going toward small scales makes it 
necessary to have a better understanding of scalar turbu- 
lence universality. To investigate the degree and the type 
of dependence of scalar turbulence on the injection mech- 
anisms we have performed two scries of numerical exper- 
iments. The advecting velocity field in both series is the 



two-dimensional NS velocity discussed in Section [II A. 
The major difference is in the injection mechanisms. As 
discussed in Section ^, the first of them is isotropic (ran- 
dom forcing) and the second is anisotropic (by a mean 
gradients g). The comparison between the two sets of 
results will permit a quantitative analysis of universality. 



A. Anisotropy 

For the anisotropic injection by a mean gradient, the 
scalar structure functions S'„(r) depend on both r = \r\ 
and 0, the angle between r and g. The latter dependence 
is best dealt with by expanding Sn{T, 4>) on a proper set 
of orthonormal functions U^''\(p) : 



Snir, 



J2slP{r)U^Hc^). 



(9) 



The most natural basis is made of the set cos{l4>), 
sin{l (/)). This is the two-dimensional version of the SO(d) 
representation group exploited in Refs. [^,|2| to analyze 
experimental and numerical turbulence data. The sine 
functions are in fact absent in the expansion (O), since 
the scalar field is statistically invariant under the trans- 
formation (j) !—>■ —(f>. Furthermore, as a result of their 
opposite symmetries with respect to the transformation 
r ^^ ~r (0 i-^ (/) + tt), even and odd orders are expanded 
as 
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(r) COS [{21 + 1) 



(10) 



(11) 



Scaling behaviors 5'„ (r) ex r'»" are expected in the in- 
ertial range. The n-th order structure function is thus 
made of a power law superposition. It is expected that 
the various contributions are ordered according to their 
anisotropic degree, i.e. Cn increases with I for a fixed n 

il- 

To extract the leading contribution to the scaling of 
5„(r), it is sufficient to project on the lowest-order com- 
ponents as 
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(12) 



(13) 



The scaling behaviors presented hereafter have been ob- 
tained by performing such angular averages over a set of 
80 snapshots equally spaced by about one turn-over time. 
We have also verified that the contributions arising from 
higher anisotropic sectors are subdominant as expected, 
i.e. that the exponents associated to 



Si^Hr)^ 
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2^ 



Sn{r,(l)) cos{l(P)d(p l>2, (14) 



are larger than those for I = or / = 1. A precise mea- 
surement of the subdominant anisotropic scaling expo- 
nents is rather delicate due to the cancellations needed 
for their statistical convergence and we shall not dwell on 
this. 

A remark on the projection procedure is in order. In 
numerical simulations, it is easy to perform angular av- 
erages such as (|l2|) or ( p^ ) . This is not generally the case 
for laboratory experiments where measuring the data for 
all possible directions (j) niight be problematic. It is then 
useful to note that the results obtained by using the pro- 
jections (111) or dll) essentially coincide with those where 
only the first subleading contributions are eliminated. 



This filtering is much more economic as it can be per- 
formed by a single measurement along (j) = 7r/4 for even 
orders and = 7r/6 for odd-order structure functions. 




FIG. 3. The 3rd and the 5th-order parallel structure func- 
tions for the injection by a mean gradient. In the inset, local 
scaling exponents dSn{r)/d\ogr. The measured exponents 
are Cs = 1-25 ± 0.04, Cs = 1-38 ± 0.07, with error bars esti- 
mated from rms fluctuations of local scaling exponents. 

Let us now present the results concerning the 
anisotropy decay rate, measured by odd-order structure 
functions. In Fig. 0, we show the 3-rd and the 5-th or- 
der structure functions. Accordingly, the skewness and 
hyperskewness coefficients of scalar differences scale as 
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the second-order scaling exponent being compatible, 
within the error bars, with the dimensional value ~ 2/3 
(see Fig. ^. Both behaviors ( p^ ) violate the dimensional 
prediction (R). Furthermore, whilst the skewness is de- 
caying, even though much more slowly than the expected 
r^/'^, the hyperskewness grows at small scales. The ef- 
fect of memory of the large-scale injection conditions 
observed in laboratory experiments is thus dramatically 
present also in our numerical simulations. The fact that 
clean scaling behaviors and exponents are measured here 
gives a very strong indication in favor of the fact that we 
are not dealing with a finite Peclet number effect. The 
most natural conclusion stemming from experiments, the 
previous and the present numerical simulations is that 
isotropy is not restored in the full classical sense, no mat- 
ter how large is the Peclet number. 

It is worth discussing in some more detail the rela- 
tion between this last conclusion and the observation 
previously made on the subdominance of higher-order 
anisotropic contributions to structure functions. The 
problem comes of course from the presence of intermit- 
tency : contrary to the classical theory, the pdf of scalar 
differences V{drT) for various separations r cannot be 
collapsed one onto another by a simple rescaling. Their 
change of shape with the separation r makes it crucial 
to specify how the anisotropy degree is quantified and 
the amplitude of the fluctuations sampled in the pdf's 



at various r's. The ratio S^/ Sj is for example decreas- 
ing with r, contrary to the hyperskewness in (pa). The 
blow-up of the latter is thus a joined effect of anisotropy 
and intermittency. The point to be remarked is that 
odd-order structure functions are all scaling with ex- 
ponents higher than those of structure functions of the 
same order but calculated with the absolute values, e.g. 
/ \5rT\^ V{5rT) d6rT . In other words, when fluctuations 
of similar amplitude are compared, anisotropic contri- 
butions are subdominant with respect to the isotropic 
ones. The same point is stressed in Fig. |[ where it is 
shown that the ratio between the anti-symmetric and 
the symmetric part of scalar increment pdf's decreases 
with r. Notwithstanding the conclusion on the absence 
of isotropy restoration in the classical sense, it is thus 
important to realize that anisotropic fluctuations become 
more and more subdominant as their anisotropy degree 
increases. 




FIG. 4. Ratio between the anti-symmetric and the sym- 
metric part of scalar increment pdf's 'P{SrT) for four separa- 
tions in the inertial range. 



B. Universality 

What is the degree of universality of scalar turbu- 
lence with respect to the injection mechanisms? Is 
the persistence of anisotropy just discussed signaUing 
that the scalar statistics is totally imprinted by the 
non-universal large-scale injection conditions even at the 
smallest scales? To answer these questions it is conve- 
nient to compare observables, such as even-order struc- 
ture functions or pdf's themselves, that are meaningful 
and non trivial for both types of injection. 

Let us first start from the scalar field pdf. Even though 
this is not a small-scale quantity, it is interesting that 
the two curves are different and the one for the mean- 
gradient injection is subgaussian, see Fig. ^. Note that 
this is not in contradiction with the exponential behav- 
ior predicted in |^,^,Q . The latter work considers in- 
deed the situation where the velocity correlation length is 
much smaller than the typical length of variation of the 
large-scale mean scalar profile. In our case, these two 
lengths are of the order of the velocity friction length 



in the structure function dominant expression Sn{r) 



scale, defined in Section III, and of the box size, respec- 
tively. They are thus comparable^ Similar subgaussian (j e"/^e~"'/^7-"/3(i/r)"/3-Cr. 
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FIG. 5. The Pdf 's of the scalar field normalized by its stan- 
dard deviation for the two cases of random forcing (upper 
curve, multiplied by a factor 10^ for plotting purposes) and 
injection by a mean gradient (lower curve). Note the pres- 
ence of subgaussian tails for the latter (Gaussian densities 
are shown for comparison as dotted lines). 

Let us move then to the small-scale statistics by consid- 
ering scalar increments. In Fig. 0, we present the proba- 
bility density functions for the two injection mechanisms : 
it is clear that they do not have the same shape and 
the same conclusion holds if we take only the symmet- 
ric part of the pdf's. A particular consequence is that 
structure functions cannot have a universal expression, 
as postulated in the classical theory. It is however di- 
rectly checked from Fig. that the scaling exponents are 
the same for the two types of injection. This implies that 
the pdf's, although having different shapes, rescale with 
r in the same way. 



constants C„. The resulting picture of universality is 
the following : structure function scaling exponents are 
independent of the injection details and thus universal ; 
scalar and scalar increment pdf's, and non-dimensional 
constants in the structure functions are not. A particu- 
lar instance is the one of odd-order structure functions, 
where the coefficients C2n+i vanish or not, depending on 
the symmetries of the injection mechanisms. This form 
of universality is weaker than the one in the original clas- 
sical KOC theory and agrees with the remark made by 
L.D. Landau about the universality in Kolmogorov 1941 
theory for fully developed turbulence. Remark finally 
that this universality framework is the same as the one 
arising from the zero mechanism in the Kraichnan passive 
scalar model [^-08i. 
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FIG. 7. Low-order even structure functions. Local scaling 
exponents are shown in the inset. The measured exponents 
are C2 = 0.66 ± 0.03, C4 = 0.95 ± 0.04 and Ce = LH ± 0.04. 



V. SATURATION OF INTERMITTENCY 
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FIG. 6. Pdf's of scalar increments normalized by their stan- 



dard deviations for three separations r 
lO"'^ in the inertial range. 
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The scaling exponents being universal, the depen- 
dence on the large-scale injection conditions shows up 



The cliffs observed in scalar fields are strikingly sugges- 
tive of quasi-discontinuities. When smaller and smaller 
molecular diffusivities are considered, the minimal width 
of the fronts shrinks with the dissipation scale, with their 
maximum amplitude remaining comparable to the scalar 
rms value. Simple phenomenology suggests that the pres- 
ence of such structures, corresponding to a local Holder 
exponent equal to zero, might induce a vanishing slope in 
the structure function scaling exponent curve. The fronts 
being the strongest possible events, this behavior should 
take place for large enough orders, whence the possible 
saturation ^„ — > const for high n's. 

The first and most natural way to investigate the satu- 
ration is to look directly at the scalar structure function 
scaling exponents. An alternative procedure consists in 
looking, for a fixed SrT, at how the pdf V{6rT) of the 
scalar differences varies with the separation r. The lat- 
ter is statistically more reliable than the former as less 



sensitive to the extreme tails of the pdf. 
equivalent to the pdf taking the form 



V{5rT)^r'^-^Q 



SrT 



1 



Saturation is 



(16) 



,1/2 



for values oi SrT > Trms = {{T^) - {Tfy'\ The tails of 
the (non-universal) Q function are bounded by the single- 
point scalar pdf, as it follows from the simple inequality 



(see, e.g., Ref. |39 



Prob{\5rT\ >t}< 2Prob{\T\ > r/2}. 



(17) 



Both procedures have been used for the case of the 
Navier-Stokes advecting velocity field, where we could 
efficiently use spectral numerical methods. They present 
the disadvantage of being sensitive to finite-size effects 
but the advantage of giving access to the whole field and 
thus to the pdf's. For the Kraichnan model, spectral 
methods are not very efficient and it is more convenient to 
use Lagrangian techniques (see Ref. [£l|,|0| for a detailed 
discussion) . The latter permit very precise measurements 
of the scaling exponents at the price of not giving access 
to the pdf's. 



A. Saturation w^ith the Navier— Stokes advecting flow 

The behaviors of the structure functions 5*10, S'12 and 
S'i4 in the presence of a large-scale scalar gradient are 
shown in Fig. || The scaling exponents („ vs the order 
n (up to the order 14) are reported in Fig. ^. The curve 
is compatible with a saturation exponent (^00 ^ 1-4. The 
convergence of the moments was inspected by the usual 
test of checking whether {5rTY^V{5rT) decays before the 
pdf V becomes noisy. The bulk contribution to the mo- 
ments 10, 12 and 14 is shown in Fig. |lO|. As for the second 



procedure on the pdf's, following ( 
the curves r-'^-=-'P{5rT). Their col 
the saturation and gives the unknown function Q. 



n we plot in Fig. ITl| 
apse is a signature of 




FIG. 8. Scalar structure functions of orders 10, 12 and 14 
for the anisotropic injection mechanism 
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FIG. 9. Measured scaling exponent ("„ 

for the Navier-Stokes advection. Error bars are estimated 
by the rms fluctuations of local scaling exponents. 




FIG. 10. Bulk contribution to the n-th moment, 
5rT"r{SvT), for n = 10, 12, 14. Here r is at the lower end of 
the inertial range. 
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FIG. 11. The pdf's V{SrT) , for the random forced injec- 
tion, for three values of r inside the inertial range of scales, 
multiplied by the factor r~''°° . The collapse of the curves in- 
dicates the presence of saturation and also gives the unknown 
function Q{SrT/Trms) in Eq. ([|). 

In Fig. |2| we show the probabilities J^j, V vs r of 
having 6rT > XTrms for various values of A. The paral- 
lelism of the curves is again the footprint of saturation. 



Explicit evidence for the universality of Coo is provided 
in Fig. |l^, where the probabilities J^^ j, "P of having 
5rT > 2.5Trms are shown both for the isotropic and the 
non-isotropic injection mechanisms. 
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FIG. 12. Probabilities /^"^ V vs r 

of having &rT > \Trms for various A. In the inside zoom, 
such probabilities are compensated by the factor r~''=° . 
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FIG. 13. Probabilities /„ ^^ V vs r of having 
5rT > 2.5Trm3 for both the isotropic and non-isotropic in- 
jection mechanisms. The parallelism of the curves reflects 
the universality of the saturation exponent ^cx) ■ 

Both procedures give therefore a strong evidence in fa- 
vor of saturation for the Navier-Stokes advecting velocity 
field. 



B. The case of the Kraichnan model 

The advecting velocity field in the Kraichnan model 
|19| is assumed to be Gaussian, homogeneous, isotropic, 
of zero mean and correlation function 



{{v{r,t) - v{0,t))^{v{r,t) - v{0,t))fj) 



Dr^ 



{d + ^-l)d^f3-^ 



r^rp 



(18) 
(19) 



where d is the dimension of space. The assumption 
of (5-correlation in time is of course far from the re- 



ality, but it has the remarkable feature of leading to 
closed equations for equal-time correlation functions 
Cn — _ (^^^ 1 ,t) ■ ■ ■ 0(rn,t)) of any order n (see, e.g., 
Refs. [[24[|4l[ ). This has permitted to rely the emergence 
of anomalous scaling and intermittency to the so-called 
zero modes, i.e. functions annihilated by the inertial op- 
erators appearing in the aforementioned closed equations 
|p6|-p8[. Anomalous zero modes can be calculated non- 
perturbatively for the case of a passively advected niag- 
netic field |4^ and for shell models a la Kraichnan [Q . 
For the case of the passive scalar, the expression for the 
second order correlation function C2 is known exactly |2J] 
and the inertial-range scaling behavior for the second- 
order structure function reads: 



S2ir;L) 



2x(0) 



C2id~l)dD 



r-Ca 



C2 = 2 - e, 



(20) 



where x(0), defined in (0), is the energy injection rate. 
The exponent ^2 coincides with the predictions based on 
dimensional arguments and ^ = 4/3 is the velocity ex- 
ponent corresponding to the KOC scaling of the scalar 
field. Exact solutions are not available for higher orders 
and analytical predictions are perturbative in three dif- 
ferent limits of the model : rough velocity fields ^ <C 1 
(Refs. ||3 g , p^ , p5[ ) ; large space dimensionalities d :$> I 
(Ref. [|37| ) and almost smooth velocity fields ^ ^ 2 
(Ref. P^). In all of them, the order n of the struc- 
ture functions enters into the correction to the dimen- 
sional scaling together with the small parameter of the 
expansion. This means that, for any non- vanishing value 
of the expansion parameter, the predictions for C„ will 
break down at large enough n and the very strong events 
associated to the extreme tails of the pdf's cannot be 
captured by such procedures. To deal with the rare fluc- 
tuations responsible for the asymptotic behavior of the 
exponents C„ at large n's, instantonic approaches have 
been borrowed from field theory ^,^. In the limit of 
very high d's, the problem could be solved Es], with sat- 
uration taking place for n > n* (the latter diverging as 
d — > 00). In the three-dimensional case, saturation has 
been phenomenologically suggested in Ref. ^] and in- 
ferred from an instantonic bound in Ref. pQ]. 

Since numerical simulations are a fortiori limited to the 
measure of finite-order structure functions, it is impor- 
tant to choose the most convenient conditions to observe 
saturation. In particular, we expect that, for a fixed di- 
mensionality of space, the orders where saturation takes 
place reduce with C2 = 2 — ^. This is due to the mecha- 
nism of cliff steepening. The latter is clearly favored by 
the presence of large-scale velocity components, inducing 
a coherent shear effect, and obstacled by small-scale ve- 
locity components, leading to effective diffusion effects. 
Since the amplitude of the former increases as 2 — ^ re- 
duces, we expect that saturation is favored by choosing ^ 
close to two. Approaching too closely the Batchelor limit 
^ = 2 is however problematic from a numerical point of 
view, due to the presence of strong nonlocal effects. The 
chosen trade-off is to consider the three cases ^ = 1.875, 



^ — 1.84, ^ = 1.75, in 3D. The exponents being universal, 
it is more convenient to consider the isotropic injection 
mechanism (||) and we shall look at the fourth and the 
sixth-order structure functions. 
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FIG. 14. The measured scaling exponent ("„'s joined by 
straight-lines, (a) i = 1.875; (b) C = 1.84 and (c) ^ = 1.75. 

A detailed description of the Lagrangian method can 
be found in Refs. [pllpO| and will not be reported here. 
We only remark the fact that it naturally allows to mea- 
sure the scaling of the structure functions S^nij", L) vs the 
integral scale L of the forcing. Physically, this means that 
the injection rate of the passive scalar variance (which 
equals its dissipation rate) and the separation r are kept 
fixed, while the integral scale L is varied. Anomalies, 
that is deviations from dimensional scaling exponents, 
are therefore measured directly (see (g)) through the scal- 
ing dependence on L of the structure functions. 

The measured C„'s are shown in Fig. |lj for: (a) 
C = 1.875, (b) ^ = 1.84 and (c) f = 1.75. (The error 
bars are estimated by the rms fluctuations of local scal- 
ing exponents over intervals of constant length logj^Q 5 in 
log-space.) The corresponding behaviors of 5*6 vs the in- 
tegral scale L are shown in Fig. Il^ for ^ = 1.875, S, = 1.84 
and £, = 1.75, respectively. The structure function 5*4 for 
f — 1.875 is shown in Fig. |6|. Similar high quality scaling 
laws have been obtained for ^ = 1.84 and ^ = 1.75. Some 
remarks about Fig. ^ are worth. It follows from Holder 
inequalities that the ^n curve for n > 6 must lie below the 
straight line joining C4 and Cg. It seems unlikely that C„ 
can be a decreasing function of n. Indeed, an adaptation 
of the argument presented in Ref . [2^ , indicates that a 
decreasing C„ would entail arbitrarily large temperature 
differences at very small scales, something unphysical, 
given the maximum principle for the advection diffusion 
equation. Note that stricto sensu this argument applies 
only to the unforced equation, but the scalar scaling ex- 
ponents are expected to be the same in the forced and in 
the unforced case (in its quasi-stationary phase). Hence, 
curve (a) in Fig. 04 and Holder inequalities strongly con- 
strain the behavior at higher n's. 
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FIG. 15. 3-d sixth order structure function Se vs L. Top: 
^ = 1.875. Parameters: r = 0.027, k = 5 x 10'^^ and number 
of realizations ^ 21 x 10®. The best-fit (solid) line has the 
slope p6 = 3C2 - 0.19 = 0.185 (±0.01). Middle: ^ = 1-84. 
Parameters: r = 0.54, k = 5 x 10^^^ and number of real- 
izations ~ 23 X lO''. The best- fit (solid) line has the slope 
P6 = 3C2 - 0.27 = 0.21 (±0.01). Bottom: ^ = 1.75. Pa- 
rameters: separation r = 0.032, diffusivity k = 1 x 10"^ and 
number of realizations '-^ 27 x 10®. The best-fit (solid) line 
has the slope pa = 3(2 - 0.42 = 0.33 (±0.015). 
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FIG. 16. 3-d fourth order structure function S4 vs L for 
5 = 1.875. Parameters: r = 0.027, «: = 5 x 10"^^ and number 
of realizations ~ 5 x 10®. The best- fit (solid) line has the slope 
P4 = 2C2 - 0.064 = 0.186 (±0.002). 

Of course, the presence of error bars makes it impos- 
sible, with a finite number of measured Cn's, to state 
rigorously that the curve (a) tends to a constant. It is 
nevertheless quite clear that the latter provides strong 
evidence in favor of saturation of the scaling exponents 
of the structure functions. From curves (b) and (c) we 
can also notice that the dependence on 1^ of the order 
n where saturation takes place is indeed increasing with 
C2 = 2 — ^, in qualitative agreement with our previous 
simple physical arguments and the result in Rcf. Q . 



VI. FRONTS 

The presence of fronts in the scalar field (See Fig. [l7| ) 
is a major characteristic of scalar turbulence. How is it 
related to intermittency saturation ? As we have shown 
in Fig. Q3, the probability to have a fluctuation which is 
0(TrmsTacross a separation r goes as a power law r''°° . 
This result admits two different interpretations which are 
associated to very different physical pictures. 

As for the first scenario, assume that the most im- 
portant contribution to excursion of 0{Trms) comes 
from jumps of the scalar field occurring across a very 
small lengthscale (comparable to the diffusive scale 77 ~ 
{k^ /e^Y^'^). We shall call these jumps mature fronts, 
to indicate that the process of gradient steepening is 
brought to its utter development. Then, the probabil- 
ity of observing large scalar fluctuations at scales larger 
than r is determined by geometrical considerations. Sup- 
pose, for instance, that the fronts have an extension 0{L) 
in the direction transverse to the jump, that is the gra- 
dients 0(Trms/ri) ^-re occurring on a set of points which 
form a line. The probability of intercepting a front across 
a distance r is then proportional to r itself. In this case, 
the saturation exponent would be equal to unity. Since 
the measured exponent for the probability is (^00 — 1-4 



we infer that mature fronts live on a set of dimension 
Dp ~ 0.6. A dimension smaller than one means that the 
largest gradients occupy a fraction of space smaller than 
that occupied by a smooth line, thus mature fronts ap- 
pear as a collection of short, broken lines where a steep 
gradient takes place (see Fig. |l8|). What is peculiar of 
this scenario is that, as the Peclet number goes to in- 
finity, the thickness of the fronts decreases, and larger 
gradients take place, but the probability of finding a ma- 
ture front remains constant. This physical picture has 
some analogies with that encountered in Burgers turbu- 
lence : there, it is well-known |^2[ that the presence of 
isolated shocks (i.e. mature cliffs of unit co-dimension), 
connected by smooth ramps, is the cause of intermit- 
tency saturation. In the context of this purely geometri- 
cal interpretation the time variability does not play any 
role, which amounts to say that mature fronts should be 
long-living objects. This remark leads us to consider a 
different possibility. 
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FIG. 17. A snapshot of the scalar field, T, for the injec- 
tion by a mean gradient directed from right to left. Colors 
are coded according to the intensity of the field: white corre- 
sponds to large positive values, black to large negative values. 

The second scenario is based on the idea that the pro- 
cess of front formation starts from a 0{Trms) fluctuation 
of the scalar field generated at large scales by injection 
mechanisms. This profile gets steepened by the action of 
large-scale converging components of the velocity field, 
until the stretching exerted by small-scale velocity fluc- 
tuations overcomes the tendency of the front to steepen. 
As a result, the front survives only down to a scale r, 
larger than the dissipative scale 77. The probability that 
the steepening process stops at a scale r decreases with 
r according to r^°' . At variance with the former picture, 
here the probability of observing a mature front vanishes 
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thus it is strictly zero for an infinite Peclet num- 



as r]^ 
ber. 

It has to be remarked that in both cases, the dimension 
of the set hosting mature fronts has to be Dp — 2 — 
^oo , but the number of points which belong to this set is 
vanishing with increasing Peclet. Thus, the observation 
of a dimension Dp — 0.6 (see Fig. |l^ ) does not allow us 
to discriminate between the first and the second picture. 
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FIG. 18. The union of the sets where mature fronts are 
located. The snapshots over which the union is performed 
are 80, regularly spaced by one eddy turnover time. 



measure the largest excursion A^j = max/^ jf^raTj occur- 
ring on a sub-interval I2 C /i of length r2. If the mature 
fronts dominate the statistics, as in the first scenario, 
a one-dimensional section of the scalar field then would 
look like a sequence of steps. The distribution of A^j 
should then be sharply peaked around Sr-^T for any r2. 
On the contrary, if mature fronts have a relatively little 
probability with respect to non-mature ones, the typical 
value for A^j moves to smaller and smaller values while 
decreasing r2 . In Figure ^ we show that the second pos- 
sibility is realized, and thus non-mature fronts dominate 
the statistics of the large excursions. 
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FIG. 20. The pdf of Arj - maximal excursion across a 
length r2 - under the condition that an excursion larger than 
2.5Trms takes place at a scale ri, for three different r2/ri. All 
distances belong to the inertial range. 
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FIG. 19. The box-covering dimension of the mature cliffs 
set. The number of boxes of size r necessary to cover the set 
represented in Fig. hq decays as N{r) ~ r~^^ with Dp — 0.6. 

In a nutshell, in the first scenario intense small-scale 
excursion of scalar field are essentially generated by ma- 
ture fronts, whereas in the second one mature fronts play 
a statistically insignificant role. 

In order to quantify the weight of mature fronts in 
the build-up of strong fluctuations and thus distinguish 
among the two options listed above, we computed the fol- 
lowing conditional probability. Given a jump of height 
i^riTj = 0{Trms) across an interval /i of length ri, we 



Summarizing, saturation comes from a self-similar pro- 
cess and the signature of the exponent (^ac is felt through- 
out the inertial range down to the dissipative scales. In 
this sense, both mature and non-mature structures carry 
the imprinting of saturation: indeed the dimensions of 
the set of mature and of the set of non-mature fronts are 
equal to 2 — Coo- Nevertheless, since the appearance of 
a mature front is a relatively rare event, the dominant 
contribution to structure functions is due to non-mature 
chffs. 



VII. PLATEAUX 

Fronts separate regions of space where a very efficient 
mixing takes place, and fluctuations of the scalar fleld are 
small. To investigate the statistical signature of these 
plateaux we computed the low-order moments of scalar 
increments, shown in Fig. ^. We observe that the scal- 
ing exponents of negative order display an almost linear 
dependence on the order, close to p/2. Thus, the scalar 
fleld within the plateaux is not smooth but it is charac- 
terized by a Holder exponent 1/2. This would entail a 
single rescaling factor, r^^^, for the inner core of the pdf 
of scalar increments, which seems to be the case, as in 
Fig. E2. This result suggests that the nature of scalar 
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turbulence within the plateaux could be understood in 
simple terms, although, to our knowledge, no definite ex- 
planation of these observations is available up to now. 
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FIG. 21. Scaling exponents of structure functions for low, 
fractional orders, in the random-forced case. Moments are 
taken with absolute values. The dotted line is p/2. 
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FIG. 22. Probability density functions of scalar increments 
for three separations within the inertial range, rescaled by 

.1/2 



VIII. SCALAR DISSIPATION 

Let us focus our attention on the passive scalar equa- 
tion ([l|), in the presence of a random forcing (0). Denot- 
ing SrT = T{x)—T{x') the scalar increment over the sep- 
aration r = a; — a;', it is easy to obtain p9| the equation 
for the even-integer structure functions S2n (odd-integer 
structure functions are trivially vanishing for isotropic 
injection mechanisms): 



dS2n{r) 
dt 



where 



and 



+ V • {W{5rTf^) = J2n{r) (21) 



2n{{SrTf''-^H), (22) 



H EE {K,{A + A')6rT\SrT) 



(23) 



(24) 



are the ensemble average of velocity increment and dissi- 
pation conditioned on the scalar increment drT, respec- 
tively. The Laplace operators in (e3) are evaluated at x 
and x' respectively. In our DNS simulations we replaced 
the Laplace operator in Eq. (m with a bi-Laplacian dissi- 
pative term. The definition of the conditional dissipation 
H is thus: H{SrT) = -k{{A^ + {A')^) SrT\6rT}. The 
choice of a laplacian or a bilaplacian is immaterial to all 
the forthcoming arguments. In isotropic conditions the 
vector W has the form W = Wr, and is even under a 
sign reversal of the scalar field, whereas H is odd. The 
information carried by the set of equations (pl|), which 
hold for any n, can be expressed in a single equivalent 
equation for the pdf of scalar differences 



dtV{5rT) + V • {WV{5rT)) = - 



dSrT 



{HV{5rT)) 



(25) 



In the statistically stationary state this equation provides 
a dynamic link among the three quantities W,H,V. Let 
us give a few examples which will be relevant to our sub- 
ject: (j), if V is exponential and H is linear, then W is 
linear; (ii), if V is exponential and H is constant, then 
then W is constant; {Hi), if P is Gaussian and W is 
constant, then H is inversely proportional to the scalar 
fluctuation. 




8r T/Trms 
FIG. 23. The longitudinal component of the W function vs 
SrT for two different values of r in the inertial range. 

Let us now look in detail at the numerical results: W 
is shown in Fig. E3. The negative part for large values of 
6rT is physically intuitive : the formation of strong cliffs 
amounts to the build-up of intense scalar gradients, and 
takes place whenever the gradient alignes with the eigen- 
vector of the strain with a negative eigenvalue. In other 
words, to obtain a large fluctuation it is necessary that 
the flow favors the approaches of elements of fluid that 
were initially far apart and characterized by very different 
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values of the scalar field. Conversely, the positive peak 
for small scalar fluctuations shows that the "plateaux" 
are regions of extremely efficient mixing, dominated by 
small-scale disordered velocity fluctuations with high tur- 
bulent diffusivity. 

The conditional dissipation function H{6rT) is shown in 
Fig. ^. The curve closely follows the linear ansatz jlj] 
up to values comparable to Tmis- For larger values the 
curve is however distorted and has a definite tendency to 
flatten out. 
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FIG. 24. The conditional dissipation H vs 5rT for two dif- 
ferent separations r in tiie inertia! range. 

As discussed in Sec. M, the tails of scalar incre- 
ment pdf's V{8rT) for various separations r in the in- 
ertial range all collapse onto a single curve : V{SrT) = 
r(^Q[A:IL), The curve r-'^'^P{SrT) is shown in Fig. H 
for one separation in the inertial range. Its tails (and 
those of Q) decay roughly as an exponential. 
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FIG. 25. The scalar increment pdf for a separation within 

the inertial range. 

The one-point pdf was shown in Fig. ra and clearly de- 
cays as a Gaussian in the random-forced case. As a con- 
sequence, the exponential tails observed in Fig. |2^ cannot 
continue for ever and a shar per than Gaussian fall off (see 
the arguments of Sec. VIII) is expected to set in. 

From the comparison of the graphs of W,H and V 
it appears that for fluctuations small compared to Trms 



we are under the conditions (i), with an almost linear 
_ff, an approximately exponential behavior for V and a 
W ~ \SrT\. This regime is characteristic of "plateaux" 
where efficient mixing occurs {W > 0) and dissipation 
is proportional to the intensity of the fluctuation. For 
larger values SrT of the order of a few Trms we can recog- 
nize the case (ii), characterized by the presence of fronts. 
Here, the dissipation is not anymore proportional to the 
scalar excursion, rather tending to be independent of it. 
It has to be remarked that, although large gradients may 
develop within a front, they have a negligible statisti- 
cal weight, so that the main contribution to the dissipa- 
tion comes from relatively small fluctuations of the order 
K^/^eg ev . The regime {iii), where the strength of 
fluctuations is dominated by the one-point pdf, is out of 
reach for our statistics: indeed, a rough estimate of the 
amplitudes where this regime sets in is 10 Trms- 

When the velocity fleld is rapidly changing in time, the 
expression for W can be derived analytically, a fact which 
allows for further considerations which are presented in 
Appendix H. 



IX. SCALAR FLUXES 

The classical theory of scalar turbulence presented in 
Section || is based on the picture of a cascade of scalar 
fluctuations from large to small scales which takes place 
at a scale-independent rate eg. More quantitatively, the 
constancy of the scalar flux is expressed by means of an 
exact relation, the Yaglom's law, which in two dimensions 
reads 



{[Srv ■ r] [SrTf) = -2ee7 



(26) 



for any r within the inertial range of scales. In the 
anisotropic case this relation still holds for the projec- 
tions onto the isotropic sector. In Fig. 26 we show the 
velocity-scalar-scalar correlation which appears in (|26| ) 
computed numerically for both the mean gradient and 
the random forced case. 
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FIG. 26. The Yaglom's law for the shear-driven (x) and 
the random-forced (-I-) case. 
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The statistics of the scalar flux displays as well a re- 
markable interniittency. To quantify it, let us intro- 
duce the fluctuating scalar flux through the scale r as 
Ilr — 6rv{SrT)'^. We consider moments dllrl'') ~ r""" , 
which have been studied in the context of experimental 
turbulence in Ref . |^^ . As shown in Figure ^, the com- 
puted scaling exponents ap deviate significantly from the 
dimensional prediction ap = p. The behavior of the ex- 
ponents for large orde rs can be explained with the aid of 
the results of Section VIII. We observed that the aver- 



age velocity fluctuation becomes nearly constant for large 
enough scalar fluctuations (see Fig. 23). In other words, 
at small scales we observe a weak correlation between 
the scalar and the velocity fluctuations. A very rough 
estimate, based on the assumption that velocity-scalar 
correlations are small, gives us ap ~ p/3 + C,2p for large 
p, which is remarkably close to the observed values. This 
approximation suggests that there will be no saturation 



for ap^ and that asymptotically Up 
experimental data reported in Ref. p 
ture. 
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FIG. 27. Scaling exponents Qp for various moments of the 
scalar flux. The dotted line is the dimensional prediction 
ap = p and the solid line is the asymptote Qp = p/2> + Coo- 



X. CONCLUSIONS AND PERSPECTIVES 

There are several turbulent systems other than passive 
scalars which are characterized by the presence of struc- 
tures : "plumes" in turbulent convection, "vortex sheets" 
in shear-driven hydrodynamic turbulence and "current 
sheets" in magnetohydrodynamics. It is thus natural to 
ask whether some or all of the results obtained for pas- 
sive scalars apply to these other systems. 
Small-scale anisotropy. A phenomenology similar to 
the one for the passive scalar has been identified for 
Navier-Stokes turbulence in the presence of a homoge- 
neous shear flow pM and for kinematic magnetohydro- 
dynamics both in the presence of incompressible fSq ] 
and compressible p6| velocity fields. Concerning Navier- 
Stokes turbulence, an anisotropic component of the vor- 
ticity third-order correlation is found to be independent 
of the Reynolds number in the presence of a large-scale 



shear [Q . A mecahnism in terms of coherent structures 
has also been proposed for the emergence of small scale 
anisotropy |5J]. Even if the effect is less pronounced 
than in the passive scalar case, it is likely that large- 
scale anisotropics do not decay down to the small scales. 
The same happens in magnetohydrodynamics turbulence 
where it has been found that for incompressible velocity 
fields the hyperskewness of the magnetic field diverges 
at small scales Q. No systematic study has been per- 
formed in the context of turbulent convection. 
Universality. There is good agreement about the univer- 
sality, with respect to energy pumping details, of scaling 
exponents in homogeneous and isotropic hydrodynamic 
turbulence. For anisotropic turbulence it appears that 
the exponents in the isotropic sector coincide with those 
found in isotropic turbulence 133]. A similar conclusion 
has been found analytically for simple models of kine- 
matic magnetohydrodynamics |pq,pq|. No results are up 
to now available for the full magnetohydrodynamics. As 
for turbulent convection, there the system is self sus- 
tained and does not need external forcing, thus the prob- 
lem has to be rephrased in terms of universality with 
respect to boundary conditions. We are not aware of any 
extensive analysis of this issue. 

Saturation. The intrinsic instability of vortex sheets 
seems to rule out the possibility of a saturation of in- 
termittency for Navier-Stokes turbulence. Conversely, 
the extreme stability and strength of plumes in convec- 
tive turbulence suggests that saturation may occur there. 
Indeed, buoyancy forces may collaborate with compres- 
sional contributions of large scale velocity leading to 
sharper and more long-lived fronts than those observed 
in passive scalar turbulence. The situation with the mag- 
netohydrodynamics turbulence remains to be clarified. 
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APPENDIX A: DOMINANT AND 

SUBDOMINANT TERMS IN THE STRUCTURE 

FUNCTIONS 

In order to have a better insight on the angular depen- 
dence of scalar structure functions, it is useful to consider 
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a simple example where the problem can be tackled an- 
alytically. This is the case for the second order structure 
function in the Kraichnan advection model |2^,|l^. In 
this model the velocity field v — {va, a = 1, . . . , d} ad- 
vecting the scalar is incompressible, isotropic, Gaussian, 
white-noise in time; it has homogeneous increments with 
power-law spatial correlations and a scaling exponent ^ 
in the range < ^ < 2 : 



(K(r,t)-w„(0,0)][w/3(r-,t) 



Mo,m 



where, 



Do.0{r) = Dor^ i^ + d-l)S^p-^ 



far 13 



(Al) 



(A2) 



The specific form inside the squared brackets follows from 
the incompressibility condition on v. 

because of the white-in- 



As remarked in Sec. VB 



time velocity field, this model leads to closed equations 
for single-time multiple-space correlation functions such 
as C„ = (6{ri,t)...9{rn,t)) at any order n (see, e.g., 
Refs. [|4[|4l)). 

Their general form in the stationary state (i.e. dtCn = 
0) maintained by the large scale gradient g is Uw : 






rj) + K6^p]d^^d^^{9{r,)...9{r,,)) 



^9agpGap{r, - rj){e{ri). ._. ._. .9{rn)) 

■ / ■ i 7 






■r,)dt^{0{r^)...,..e{r,,)), (A3) 



with Gapir) = {va{r,t)vi3{0,t)) and where the transla- 
tional invariancc of scalar correlation functions has been 
exploited. For the sake of simplicity, let us focus the 
attention on the equation for the two-point correlation 
function C2 (r, (f>) in the two-dimensional case, with (j) de- 
noting the angle between the mean gradient g and the 
separation r. Restricting ourselves to the inertial range 
of scales, such equation has the simple form 



9 



r 

-5V[(C- 



e + 1 



di 



C2{r,(t>) 



I) 



Ccos 



(A4) 



Among the possible functional basis through which 
C2{r,4>) can be decomposed, the S0(2) representation 
functions |Q t/^'^ (6) = exp(— «/</)) turns out to be par- 
ticularly useful pq]. Indeed, exploiting this representa- 
tion we have 



C2(r,0)-^cf)(r)cos[(2O0], 



(A5) 



i=0 



(only even values of / are involved due to the invariance 
of C2 {r, (j)) under the transformation ^-> + tt) and the 



equations for each of the projections Cj {r) are closed 
and independent, and can thus be separately solved. It 
is easy to show that the solution for the isotropic com- 
ponent C2 {r) is simply : 



C2 {r) = const 



, y -^ + ^r2, (A6) 

2(2-0 4 ' ^ ' 



where the dimensional constant appearing in the right 
hand side is the only non-trivial homogeneous solution 
(zero mode), of the equation for the second order correla- 
tion function. Here, L is the integral scale of the problem. 
It then follows that the structure function in the isotropic 
sector Z = is proportional to 



_2!iLr2-« 



2(2-4)' 4 

sum of power laws. Conversely, the structure function 
((r(r) - T(0))^) for the total field T{r) = Qir) + g ■ r in 
the isotropic sector goes as : 



Sfir) 



9 



2L« 



.2-? 



2(2-0 



(A7) 



This means that the two structure functions, for the fiuc- 
tuating field 9{r) and for the total field T{r) behave, in 
the isotropic sector, in the same way ex r^~^, with de- 
creasing r; but while those for 9{r) are sums of power 
laws, structure functions of T{r) have a pure power law 
scaling. It is then useful to work with the total field T{r) 
to obtain cleaner scaling behaviors. 



APPENDIX B: SCALAR DISSIPATION AND 
SATURATION IN THE KRAICHNAN MODEL 

Let us now specialize the analysis (see also Ref. ^^ ) for 
the Kraichnan advection model (Al). Exploiting homo- 
geneity, incompressibility and the (5-correlation in time, 
V • (W{5rT\^'') can be evaluated exactly and its expres- 
sion reads 119(1: 



V.(W^(5.T)-) = --A_|. 



rj(r 



dS2n{r) 
dr 



(Bl) 



where rj{r) = D{){d — l)/2r^ is the relative diffusion op- 
erator. On the contrary, J2n{f) does not involves solely 
two-point averages and this means that, in general, it 
cannot be simply expressed in term of structure func- 
tions. 

In the stationary state (i.e. dS2n/dt = 0) and for separa- 
tions, r, in the inertial range of scales (i.e. rjg <^ r <ti L) 
a power-law behavior for S2n{T) (i.e. S2n{'r) = 5271?"^^") is 
expected. The same law, with exponent C,2n — C21 holds 
for J2n to ensure the balance between the advective and 
the diffusive term, J2„, in Eq. (|2l|). Imposing the balance 
in Eq. ( pT| ) for power law coefficients and using (Bl), the 
following equation is obtained: 



C,2n{d+C,2n - C2) = dC,2—^ — , 
S2n 32 



(B2) 



15 



where we have posed J2„ = J2n?'^^"~^^- 

Let us discuss the consequences of saturation on 
Eq. ( pl| ) and thus on (B2). When scahng exponents sat- 
urate, for n >_n we have C,2n = Coo and the left hand 
side of Eq. (B2) does not depend on n any longer. To 
ensure the balance for all n >n, the same independence 
on n thus holds also in the right hand side of (B2), that 
implies; 



.?2n 52 
S2n J2 



C 



2n((5^r)2"-iiJ) ^ CJ2 



So 



(B3) 



where C may depends on ^ and/or d, but not on n. No- 
tice that complete saturation, i.e. Coo — C2, yields C — I. 
We now focus on the consequences of (83) on the behav- 
ior vs SrT of H{6rT) and K{SrT), with 



K = K{\V6rT + VSrT\ \6rT), 



(B4) 



a conditional mean related to dissipation. 

To do that, we have to assume a specific form for the pdf 

of scalar increments P{drT). Our attention being focused 



SrT 



on values oi6rT > Trms, so that V{6rT) oc r'^- Q T^ 

statistical relations involving high-order structure func- 
tions (e.g. their ratio) are thus fully controlled by Q. As 
a first choice, fully justified solely for the extreme tails 
SrT ^ Trms, wc Can thus assume a Gaussian shape for 
Q. When doing this, recalling that for a Gaussian pdf 



we have {(SrT)^") = (2r 
realize that the second of 



(SrT)"^^ ^), we immediately 
^) is satisfied by 



H = C- 



J2 



SrT S2 



for 



SrT ^ Tr; 



(B5) 



It is now possible to show that, as a consequence of ho- 
mogeneity alone (see Ref. [|o|) H and K are related by 
the following identity: 



HP 



d 



dSrT 



[KP]. 



(B6) 



from which, using the Gaussianity of Q, one obtains that 
at large values oi SrT : K oc j^. 



It immediately follows according to (B5) that the condi- 
tional mean K behaves as 



^"j|^°^(^- 



(B7) 



We have already stressed that the behaviors vs SrT of 
the conditional means H{SrT) and K{SrT) depend on 
the specific assumption made on the asymptotic shape of 
Q. However, in most experimental situations, one does 
not have access to the strongest events controlling the 
extreme tails of Q. Although one can clearly observe the 
rescaling ma) and thus gather information on the satu- 
ration, Q might still not have relaxed on the single point 
pdf. 
To better illustrate this point, let us assume that we can 



control scalar increments statistics for SrT > Trms, but 
far from the asymptotic situation SrT ^ Trms- We thus 
rule out Gaussianity for Q and assume a slower decaying, 
e.g. of exponential type, as suggested by the numerical 
results. When doing this, recalling that for an exponen- 
tial pdf we have (((JrT)^") = 2n((|(5^r|)2"-i), we imme- 
diately realize that the second equation in (B3) is now 
satisfied by 



H = sign{SrT) 



CJ2 



for 



SrT > Tr: 



(B8) 



From the identity (B6), it also follows that K is now oc H. 
We can conclude that, at least in the range SrT > Trms 
(but still far from SrT ^ Trms) and for an exponential 
decay of Q in (Cq), saturation is thus compatible with 
the behaviors of H and K displaying no dependence on 

SrT. 
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